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Abstract 

Using perturbative renormalization group we study the influence 
of random velocity field on the kinetics of single-species annihilation 
reaction A + A — > at and below its critical dimension d c = 2. We 
use the second-quantization formalism of Doi to cast the stochastic 
problem into a field-theoretic form. The reaction is analyzed near two 
dimensions by means of two-parameter expansion in e, A, where e is 
the deviation from the Kolmogorov scaling and A is the deviation from 
the space dimension 2. All of the relevant quantities are evaluated to 
the second order of the perturbation scheme and long-time asymptotic 
behaviour of stochastic system under consideration is analyzed. 



1 Introduction 

The annihilation reaction A + A — > 0, where the reactans are influenced 
by some external advective field, are one of the simplest examples of the 
nonlinear statistical systems. This type of reactions may be observed in 
diverse chemical, biological or physical systems. Pollutants in the atmosphere 
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of the Earth is typical example of such stochastic system. For this type of 
reaction in low space dimensions the usual description by means of kinetic 
rate equation is not sufficient and the effect of density fluctuations must be 
taken into the account. The origin of this effect could be tracked down to 
the physics of diffusion [TJ [2] . 

A theoretical analysis of this reaction has been carried out for the case of 
time-independent random drift [3J. However, for a more realistic description 
of the effects of the velocity fluctuations a time- dependent random drift would 
be more appropriate. The most suitable approach to this end is based on 
the use of the stochastic Navier-Stokes equations [3J. To gain an insight into 
the problem, it is useful to consider for the velocity fluctuations a simplified 
model less complicated calculationally. In this work the Kraichnan model 
with finite correlation time [5] is used for modelling the advective velocity 
field. The motivation for using it is not only its simplicity, but this model 
also reproduces many of the observed phenomena in the genuine turbulence 
mass transfer [Bj. In this model there are two special limits: the self-similar 
white- in-time velocity field ("rapid-change mode") and the time- independent 
(quenched) velocity field. Both of them are considered in this paper. 
In this paper the long-time behaviour of the annihilation reaction A + A — >■ 
is studied in the presence of a synthetic velocity field at and below its critical 
dimension d c = 2. The second-quantization formalism of Doi |7J is used 
to cast the reaction part into a field-theoretic form. It is assumed that 
the statistics of velocity field is Gaussian with finite correlation time. The 
theoretical model for such a reaction is constructed and with the use of the 
perturbative renormalization group possible IR regimes are calculated to the 
second order. 

2 Field-theoretic model 

A classical reaction system may be cast into the language of second-quantization 
with the use of the creation and annihilation operators ip and ^ and the vac- 
uum state |0) [7] satisfying the usual commutation relations, 

[V(x),V» t (x)] = tf(x-x) 

[V(x),V(x')] = [V' t (x),^(x)]=0 
^(x)|0)=0,(0|V t (x) = 0,(0|0) = l (1) 
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Let P({rii},t) represent the joint probability distribution function for ob- 
serving rii particles at positions x,, then the state vector of the classical 
many-particle system is defined as the sum over all occupation numbers 

l^(*)> = E p « n *>,*)l{M>, ( 2 ) 

m 

where the basic vectors are conventionally defined as 

iK}) = n^^)no). (3) 

i 

Now the whole set of coupled equations for the probability density functions 
can be rewritten in the compact form [H [7] 

^ t mt)) = -Him), (4) 

where H = Ha + Hd + Hr. It can be shown [8] that for the annihilation 
reaction studied these terms are expressed as 

H A = j rfx^V[v(x,t)V(x)], 

H D = -D j dx^V^x), 

H R = X D j dx(^) V, (5) 

corresponding to the advection, diffusion and reaction part. Because of di- 
mensional reasons we have extracted the diffusion constant Dq from the rate 
constant K = XqD . Mean values of physical quantities may be expressed 
j8] via the relation 

(A(t)) = (0\T(A{[^{t) + l}i>(t)}eM- I" 'fydt + rio [ dx^x, 0)) ) |0) (6) 



The expectation value of time-ordered product (EJ) may be cast by the stan- 
dard operation [9] into the form of a functional integral over complex-conjugate 
scalar fields ^(x, t) and ^(x, t): 

(A(t)) = [ V^V^A{[^(t) + l]^(t)}e Sl (7) 



where the unrenormalized action Si for the annihilation reaction A + A — >■ 
is 

Si = - J dt J dx{^d t ip + V f V(W0 - Aj^V 2 ^ + 

X D [2^ + (V^V + n Q [ dx^(x, 0)}. (8) 



The most realistic description of the velocity field v(x) is based on the use 
of stochastic Navier-Stokes equation. However, in this paper we shall study 
a simplified model in which we prescribe statistical properties of the velocity 
field. Let assume that v(x) is a random Gaussian variable with zero mean 
and the correlator 
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—^P i:j (k)D v (u, k) ex V [-iu(t - t') + ik.(x - x')]. (9) 



is 



Here -Pij(k) = 5ij — kikj/k 2 is the transverse projection operator, k = |k 
the wave number and the kernel function D v is assumed to have the following 
form 

, g D 3 k 2 - 2A ~ 2 ' , . 

D ^ - JVuiDi^r ^ 

Here go is the coupling constant (small parameter of the ordinary perturba- 
tion theory) and the exponents e,A and rj play the role of small expansion 
parameters. They could be regarded as an analog of the expansion param- 
eter e = 4 — d in the usual sense of dimensional regularization. However, in 
this paper e should be understood as deviation of exponent of the power law 
from that of the Kolmogorov scaling |6], whereas A is defined as the devia- 
tion from the space dimension two via relation d — 2 + 2A. The exponent 
rj is related to the reciprocal of the correlation time at the wave number k. 
The parameter uq may be used for labelling of the fixed points and has the 
meaning of the ratio of velocity correlation time and the scalar turnover time 
[5]. Although in the real calculations e is treated as a small parameter, the 
real problem corresponds to the value e = 4/3. 

It is interesting to note that the model for the advection field v(x) contains 
two cases of special interest: 

(a) in the limit u — > oo, g' = go/u^ = const we get the 'the rapid-change 
model' D v (u, k) — > g' Q Dok~ 2 ~ 2A ~ 2e+ri , which is characterized by the white- 
in-time nature of the velocity correlator. 
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(b) limit u Q — > 0,g f Q = go/u = const corresponds to the case of a frozen 
velocity field D v (u,k) —> gQD 2 7i5(oj)k 2A ~ 2e , when the velocity field is 
quenched (time-independent). 

The averaging procedure with respect to the velocity field v(a;) may be per- 
formed with the aid of the following action functional 

S 2 = -~J dt dt' dx d-x' v(t, x)D;\t - t', x - x')v(t', x'), (11) 

where D~ l is the inverse correlator (J9]) (in the sense of the Fourier transform). 
The expectation value of any relevant physical observable may be calculated 
using the complete weight functional W(ip \ ip, v) = e Sl+S2 , where Si and 5*2 
are the action functionals ([HD and ffTTT) . 

3 Power counting and UV renormalization 

In order to apply the minimal subtraction scheme for the evaluation of 
renormalization constants, an analysis of possible superficial divergences has 
to be performed. For the power counting in the actions ([8]) and ffTTl) we 
use the scheme [4], in which to each quantity Q two canonical dimensions 
are assigned, one with respect to the wave number dq and the other to 
the frequency dq. The normalization for these dimensions is conventional: 
d 1 ^ = —df = 1, d\ = —d k x = 1 and d% = d k = 0. The canonical (engineering) 
dimensions for fields and parameters are derived from the condition for ac- 
tion to be a dimensionless quantity. 

The quadratic part of the action ([H]) determines only the canonical dimensions 
of the product of fields ip^ip. In order to keep both terms in the nonlinear 
part 

A A) J dtdx[2^ + (^) 2 }ip 2 

of the action, the field ip* must be dimensionless. If the field tp* has a positive 
canonical dimension, then the quartic term should be discarded as irrelevant 
by power counting. The action with cubic term only, however, does not 
generate any loop integrals corresponding to density fluctuations and thus is 
uninteresting for the analysis of fluctuation effects. 

Using this choice we get the following canonical dimensions for fields and 
parameters in d- dimensional space: 
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Here, c/q = 2dg + d% is the total canonical dimension and it is determined 
from the condition that the parabolic differential operator of the diffusion 
scales uniformly under the transformation k — > fik,u — > [j?uj. 

The model is logarithmic for e = A = 0, the ultraviolet (UV) divergences 
have the form of poles in various linear combinations of e and A. The to- 
tal canonical dimension of an arbitrary one-particle irreducible Green (1PI) 
function T = (<&... $)i-i r is given by relation 

d T = d + 2 - N^, (12) 

where N$> = {Nm, N v } are the numbers of corresponding external fields. 
Superficial UV divergences may be present only in those T functions for which 
dr is a non-negative integer. Although the canonical dimension of the field ^ 
is zero, there is no proliferation of superficial divergent graphs with arbitrary 
number of external ip^ legs- This is due to the fact that n^t < n^, which is 
easily seen by a straightforward analysis of the graphs pQ. As has already 
been shown [5] that the divergences in (1PI) Green functions containing 
at least one velocity field v could be removed by the only counterterm of 
the form if)^d 2 ip, which leads to the following renormalization of parameters 
go, D and u : 

D = DZ D , g = gfi 2e+v Z g , u = u[i v Z u , (13) 

where \i is the reference mass scale in the minimal subtraction (MS) scheme. 
The renormalization constants Z g , Z u and Zo satisfy the relations 

Z g = Z~ D \ Z u = Z~ D \ (14) 

It should be noted that the graphs corresponding to r^t^ and r^t^t^ differ 
only by one external vertex and thus give rise to equal renormalization of the 
rate constant Ao-Do- Singularities presented in these (1PI) Green functions 
may be eliminated by the following renormalization of Ao: 

A = A/x- 2A Z D %. (15) 

The explicit form of the renormalization constants Zo and Z\ at the two-loop 
order can be found in the appendix. 
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4 Fixed points 



The coefficient functions of the RG operator 



n 9 



dfj, ^ dgi d D 

9i 



where the bare parameters are denoted with the subscript "0", are defined 
as 



d\nZ D 
Id = fi- 



n 9i 
, Pi = /I — 



(17) 



with the charges gi = {g,u, A}. From this definition and from relations 
(111,(111 and ([T5D it follows that 

P g = g(-2e- V + 3 lD ), f3 u = u[- V + 7d ], /3 a = A(2A - 7a + 7d ) (18) 

The scaling regimes are associated with the fixed points of the corresponding 
RG functions. The fixed points are defined as such points g*,u*, X* for which 
the (3 functions vanish 

p g {g\ u*, A*) = p u (g\ it*, A*) = p x (g\ u*, A*) = 0. (19) 

The type of the fixed point is determined by the eigenvalues of the matrix 
Q = {Qik = d(3i/dgk}, where is the full set of functions ([TBI and g^ is 
the full set of charges {g, u, A} . The IR asymptotic behavior is governed by 
the IR stable fixed points, for which all eigenvalues of Q matrix are positive. 
It is easy to see that the functions (3 g and (3 U satisfy relation f3 g /g — 3(3 u /u = 
2(77 — e). This means that they cannot be equal zero simultaneously for the 
finite values of the charges g and u. The only exception is the instance e = 77, 
which should be studied separately. For general case e ^ f) we have to set 
either u = or u = 00 and rescale g in such a way, that 773 remains finite 

ma- 

In what follows we present the results for fixed points, anomalous dimensions 
and eigenvalues of the Q matrix to the second order of perturbation theory. 
However, we would like to stress, that the form of ft functions ffl8|) allows to 
calculate the anomalous dimensions 773 and 7a exactly (without any second- 
order correction). 

In [5] independence of the renormalization constant Zd on the exponents 
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rj at the two-loop approximation has been conjectured. It implies that we 
may use the choice r] — 0, which we have applied in our calculations of the 
renormalization constants Zd and Z\. 

Let us consider the "rapid-change mode" (u —> oo). It is convenient to in- 
troduce new variables w = 1/u, g' = g/u 2 and the corresponding (3 functions 
obtain the form 

(3 g , = g '[ v -e + lD }, P w = w[ V - lD }, (3 X = A(2A- 7a + 7d ). (20) 

The "rapid-change model" corresponds to the fixed point with w* = 0. In 
this case four stable IR fixed points are realized: 

FP 1A: w* = 0, g'* = 0, A* = 

7d = 0, 7 A = (21) 

fii = rj - 2e, fi 2 =V, ^3 = 2A 
FP IB: w* = 0, g'* = 0, A* = -4ttA 

7 d = 0, 7 A = 2A (22) 

fii = rj - 2e, fi 2 =77, ^3 = -2A 
FP 2A: w* = 0, g'* = 8vr(2e - 77), A* = 

lD = 2e - 77, 7A = (23) 

fii = 2e - 77, n 2 = 2r]- 2e, fi 3 = 2A + 2e - 77 
FP 2B: u?* = 0, = 87r(2e-r7), 

A* = -2tt(2A + 2e - 77) - 7r 2A+ e 2e ~ r? (2e - 77)(A + e) 

7 d = 2e - 77, 7 A = 2A + 2e - 77 (24) 
Q 1 = 2e- v, ^2 = 27]- 2e, Q 3 = -2A - 2e + 77. 

For the analysis of the regime u — > (quenched velocity field) we introduce 
the new variable g" = g/u 2 . Hence the corresponding (3 functions have the 
form 

P gH =g"[-2e + 2 lD ], $ u = u[-r, + lD ], f3 x = A(2A - 7a + lD ) (25) 
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Also in this case there are four possible IR stable fixed points: 
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4A + 3e n 
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(26) 
(27) 

(28) 



(29) 



(30) 



where Q = 1.159 is a numerical constant resulting from a simple but cum- 
bersome two-loop integral. 

In the special case e = i] the functions (3 g and (3 U become proportional and 
this leads to the degeneration of fixed point. Instead of just plain fixed point, 
we have a line of fixed points in the g — u plane. 



FP 5A: A* = 

Id = e, 7a = (31) 
fti = Gi(u*)e + G 2 (u*)e, ft 2 = 0, ft 3 = 2A + e 

FP 5B: A* ^ 

lD = e, 7a = 2A + e (32) 
ft x = Gi(u*)e + G 2 (M*)e, ft 2 = 0, ft 3 = -2A - e 
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where the functions G\ and G 2 are given as 




u(3u 3 + 15m 2 + 26m + 18) In 



u 2 + 2u 
(1 + uf 



(33) 



The zero eigenvalue Vt 2 is connected with the existence of a marginal direction 
in the g — u plane (along the line of the fixed points). 

The "real problem" corresponds to the value e = 4/3, which leads to the 
famous Kolmogorov "five-thirds law" [6] for the spatial velocity statistics. 
By direct numerical calculation it may be easily shown that in this case 
(e = 4/3) the eigenvalue Q\ is positive. Of course this holds just for the 
physically possible values of parameters u and £. One has to consider u > 
and £ < 0, which corresponds to the space dimensions lesser than two. The 
conclusion is that all fixed points on the line are IR stable. 

5 Conclusions 

This paper is devoted to the study of reaction kinetics of the annihilation 
reaction A + A — > 0. In order to use the technique of the perturbative renor- 
malization a field-theoretic model is constructed. All the calculations were 
performed to the second order of the perturbation theory. The IR stable fixed 
points, dimensions and corresponding regions of stability of fixed points are 
calculated. 

The technically relatively simple model of velocity fluctuations used here is 
a convenient starting point for more realistic high-loop calculations based on 
the stochastic Navier-Stokes equations. These calculations are almost com- 
pleted and we hope to publish them in the near future. 
The work was supported by VEGA grant 0173 of Slovak Academy of Sciences, 
and by Centre of Excellency for Nanofluid of IEP SAS. This article was also 
created by implementation of the Cooperative phenomena and phase transi- 
tions in nanosystems with perspective utilization in nano- and biotechnology 
projects No 26220120021 and No 26220120033. Funding for the operational 
research and development program was provided by the European Regional 
Development Fund. 
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Appendix A: Explicit form of renormalization 
constants Zd and Z\ 



For brevity we use the following normalization of charges 
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The following results are presented directly for the space dimension d = 2 
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v/2(l +u) -* 2 " " ^2(1 + u) -z 2 
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